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Using the electromagnetic form factor of a quark as a working example, we de- 
scribe a subtraction technique to treat infrared sensitive regions in non-inclusive 
processes. 



Inclusive hard-scattering processes, characterized by a single large mass 
scale, are investigated in QCD by using asymptotic freedom and factorization 
theorems □. But the application of QCD to the study of multiparticle final 
states, involving several mass scales, is much subtler. The main practical tool 
is provided by Monte Carlo event generators, modeling parton shower and 
hadronization. In these event generators the theory does not yet go system- 
atically beyond the leading logarithms. To incorporate next-to-leading order 
QCD corrections in parton showers, extensions of the factorization theorems 
are necessary, for which new more precise methods are needed. 

An important step in this program is to show how to decompose Feynman 
graphs into terms associated with particular regions in loop momentum space. 
In the case of a Monte Carlo event generator simulating the exclusive structure 
of the hadronic final states, the observables being computed-are not infrared 
and collinear safe. It is important to develop techniques BHIj such that even 
for such observables the integrands to be associated with the ultraviolet region 
are integrable functions. — and can in particular be integrated numerically 
through a Monte Carlou. 

Let us consider as an example the virtual corrections to the electromagnetic 
form factor of a quark (Fig. 1). The theory for this process is well known. See 
for instance the review in Ref.LJ. To simplify the calculations while retaining 
all the ingredients that are essential for our discussion, let us work in a massive 
abelian theory with scalar quarks. We denote the quark mass by m and the 
gauge boson mass by m g . We work in a center-of-mass frame in which the 
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incoming quark momenta pa and ps are in the 
2p\pg = Q 2 . We consider the amplitude 



and — z directions, with 
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where UV indicates the MS counterterrrufor the ultraviolet divergence. 

Standard power counting arguments^ determine the.pegions of momentum 
space contributing to the leading power behavior of Tan: 1) the soft region, 
where all components of k^ are much smaller than Q: ~ AQ, with A small; 

2) the p^-collinear region: k + ~ Q, k± ~ A Q, k~ ~ A 2 Q, with < k + < 

3) the ps-collinear region: fc~ ~ Q, k±_ ~ A Q, fc+ ~ A 2 Q, with — < fc~ < 0; 

4) the hard region, where all components of k^ are of order Q. 

+ 
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Figure 1: (a) One- loop graph for the electromagnetic form factor; (b) geometry of the infrared 
sensitive regions in momentum space. 

Rcf.i constructs a decomposition of T into a sum of terms, one for each of 
these regions (Fig. 2), 



T = S + A + B + H + nonlcading power , 
2 
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satisfying the following requirements: 

i) The splitting between the terms is to be defined gauge-invariantly: we 
demand that the terms be obtained from matrix elements of gauge-invariant 
operatorsLD'El. 

ii) In particular, the necessary cut-offs on rapidity integrations should be 
gauge-invariant. As wa will see, this involves the use of Wilson lines along 
non-lightlike directions □. 

iii) The evolution equations Bu'EHl with respect to these cut-offs should be 
simple, in the sense that there should be no power-law remainder terms. All 
the power-law corrections are associated with the initial construction of the 
terms S, A, B, H in Eq. (§). 

iv) The integrand associated with the hard region should be an integrable 
function even when the physical observable being computed is not infrared safe 
in perturbation theory. 

The strategy we use to construct such a decomposition is similar to the 
R-operation techniques for renormalization. See Ref.OJ for a related approach. 
We proceed from smaller to larger regions (Fig. 2). For each region, we con- 
struct a term that, added to the terms for smaller regions, gives a good leading- 
power approximation to the original amplitude in that region, and does not 
receive leading contributions from regions that are smaller or have an overlap 
with the region being treated. _ 

To see how this works, let us look at the form of the resultH for one of the 
terms. The term S associated with the soft region is 

-*s 2 r JU+ „_ , 2l . 1 
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The first term in the square brackets is just obtained by taking the soft ap- 
proximation to Eq. ([!]). This term still has singularities from the ultraviolet 
and collinear regions. The ultraviolet singularity is to be dealt with by the 
standard subtractive renormalization procedure. We treat the collinear sin- 
gularities in a similar fashion: the next two terms in the square brackets are 
subtractions terms designed to cancel the collinear contributions. 

To define these terms we have introduced two vectors ua — {u\ , u~ A , 0) , 
ub = (ug, Ug,0), lying along directions away from the light cone. The second 
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Figure 2: Soft (S), collinear (A, B) and hard (H) contributions. 



term in the square brackets subtracts the divergence from the region collinear 
to pa, i-e., k~ /k + — > 0. The non-lightlike vector ub in this term provides 
a cut-off on the region of small k + . Similarly, the third term subtracts the 
divergence from the region collinear to ps, with the vector ua providing a cut- 
off on the region of small k~ . Further inspections of the contour integrations 
in Eq. (^) for k + k~ <C k 2 shows that ua and ub must be spacelike, u\, Ug > 0, 
u Ai u ~b < 0. Note that the collinear-to-pA subtraction term has no collinear- 
to-pB singularity; indeed it is power suppressed in this region. The same is 
true with A and B exchanged. 

The important point is that the cut-offs thus introduced are defined gauge- 
invariantly: the counterterms in Eq. (^|) can be obtained from matrix elements 
of path-ordered exponentials of the gluon field along non-lightlike lines. For a 
generic direction n, define 



V q (n) = P exp yig J dzA(zn)-n 
Vg(n) = Vexp I -ig I 

\ J — ex 



dzA(zn)-n) . (4) 



Consider the product of vacuum expectation values 

~ = (0\V q (pA) V 9 {pb)\0) (01^(^)10) {0\V 9 (u B )\0) (5) 
(0\V q (pA)Vq(u B )\0)(0\V q (u A )Vq(p B )\0) ' [) 

At one loop the first factor in the numerator gives the unsubtracted soft term 
in the first line of Eq. (||), while the two factors in the denominator, involving 



4 



Wilson lines along spacelike directions, give the collinear subtractions. The 
remaining factors in the numerator cancel factors of a complete external prop- 
agator for the Wilson line. Given the one-loop formulas, this result appears 
to be unique, if we simply assume that the quantity which we calculate is the 
product of vacuum expectation values of some Wilson line operators. 

The vectors ua, ub introduced by the subtractions are not physical pa- 
rameters. Their utility comes from the fact that evolution equations in ija = 
(1/2) In |m^/u7|, yg = (1/2) In \ug/ug\ can be applied to the terms in Eq. (||) 
to extract effects associated with large logarithms EHiiJ. One of the advantages 
of the subtraction procedure-described here is that the corresponding evolution 
equations are homogeneous El. This can be contrasted, e.g., with the case of 
Rcf. 13, where the evolution equations have power-law corrections. The sim- 
pler structure of the equations may be helpful in more complicated cases, such 
as the factorization needed for the inclusion of next-to-leading corrections in 
Monte Carlo event generators. 

The collinear terms A, B in EcL-ffl) can also be given an operator definition 
in terms of spacelike Wilson lines □. The hard term H is obtained by-taking 
the massless approximation to V — A — B — S. The result for H reads □ 



H = 



8tt 2 
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where Q 2 — 2p\p B and 



l-k 2 /Q 2 
R 



In 



1 + R 

l-R 



UV 



(6) 



R 



y/l -4k 2 /Q 2 if 4k 2 /Q 2 < 1 
i^/4k 2 /Q 2 - 1 if4k 2 /Q 2 >l 
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H is independent of the choice of the vectors ua, ub- As a result of the infrared 
subtractions, the k 2 integration in Eq. (^) is regular at small k 2 . The large k 2 
behavior is to be dealt with via an ultraviolet counterterm. 

In conclusion, the subtractive procedure that we have applied gives a finite 
coefficient for the hard part of the form factor, while the counterterms have a 
simple meaning in terms of gauge-invariant operators and obey homogeneous 
evolution equations^ Similar procedures can be defined in real emission pro- 
cesses for collinear u and soft □ contributions. Calculational schemes of this 
kind will be needed to improve the accuracy of Monte Carlo calculations for 
multiparticle final states beyond the leading order. 
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